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Abstract 

We develop a method for calculating the doublet splittings of select ground-state 
single A-hypernuclei. This hypernuclear spectroscopy is conducted by supplement- 
ing the self-consistent single-particle equations with an effective interaction, which 
follows directly from the underlying lagrangian, to simulate the residual particle- 
hole interaction. Our previous investigation, performed using only the leading-order 
contributions to the particle-hole interaction, was inadequate. In the present work, 
this method is improved upon by increasing the level of truncation in the residual 
interaction to include gradient couplings to the neutral vector meson, and thereby 
incorporating the tensor force into the calculation (which is known to play a crucial 
role in these systems). As a result, we obtain a realistic description of the effect of 
the tensor couplings on the doublet orderings and splittings. 
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1 Introduction 



In the Kohn-Sham framework, the nuclear many-body system is reduced to 
a set of single-particle equations with classical fields [1,2,3]. This framework 
allows one to reproduce the exact ground-state energy, scalar and vector den- 
sities, and chemical potential, provided that the mean-field energy functional 
is accurately calibrated. In some cases, however, the single-particle levels are 
actually weighted averages of multiple states [4]. To illustrate this point, con- 
sider the ground-state of 15P17, which has the configuration (2si/2)p(ld3/2)n- 
The angular momenta of the valence proton and neutron couple; therefore, the 
calculated Kohn-Sham ground-state is actually a doublet. To determine the 
true level orderings and splittings, one can supplement the Kohn-Sham equa- 
tions with a Tamm-Dancoff Approximation (TDA) analysis of the particle- 
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hole states. The particle-hole matrix elements are sums of two-body Dirac 
matrix elements, and the particle-hole interaction is determined by the under- 
lying effective field theory. If retardation is neglected, the interaction is given 
by (Yukawa) meson-exchange potentials (with appropriate spin-isospin oper- 
ators). This approach has been used to study ordinary nuclei [5] and single 
A-hypernuclei [6]. The case of single A-hypernuclei is particularly interesting 
as no single isovector coupling is allowed, and there are no exchange contri- 
butions, since the A and the nucleon are distinguishable. This will allow us to 
focus on isoscalar exchange in the effective interaction. 

There are no free parameters in this TDA analysis. As a result, there are 
only three possible ways to adjust this method: vary the level of truncation in 
the underlying lagrangian, introduce additional degrees of freedom, or include 
higher-order contributions in the particle-hole interaction. It is of interest to 
investigate the effect of these modifications on the accuracy of this approach, 
particularly the inclusion of higher-order terms involving gradient couplings. 
We have applied this specific improvement to calculations of ground-state, 
A-particle-nucleon-hole splittings in single A-hypernuclei, such as O. The 
higher-order terms that are of interest here are those involving the tensor 
coupling to the neutral vector meson; these terms incorporate the tensor force 
into the calculation, which is known to play a crucial role in these systems 
[7,8] and does not enter in leading-order contributions [6]. 

In recent years, a number of new avenues have opened to study hypernuclei 
with increased accuracy. Of particular interest to this work are 7-ray coinci- 
dence and (e, e'K^) experiments. Recent 7-ray experiments have measured the 
fine structure of a number of light hypernuclei [9,10], including the measure- 
ment of the ground-state particle-hole splitting in j^^O. High precision (e, e'fC+) 
experiments have measured, or are set to measure, a number of similar states 
in hght hypernuclei [11,12], including the ground-states of and jfN. Un- 
fortunately, most of the states that have been measured with these techniques 
thus far lie below the range of A accessible to the Kohn-Sham approach used 
in this work; therefore, the experimental constraints on heavier hypernuclei 
are confined to upper bounds provided by (7r+,i^+) reactions [13,14,15]. It is 
this region of medium to heavy single A-hypernuclei that this work seeks to 
investigate. 

A number of recent calculations have tackled this problem. Shell model calcu- 
lations in p-shell hypernuclei were conducted using two-body matrix elements, 
accurately describing the known data [7,8]. The influence of zero-range effec- 
tive ANN interactions on p-shell hypernuclei has also been investigated [16]. 
Another model of interest uses strangeness changing response functions to 
calculate the spectra of y O and y Ca; the resulting ground-state particle-hole 
splittings are small [17]. The spectra of O has also been calculated from a 
folded diagram method using realistic hyperon-nucleon potentials [18]. 
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In Section 2, we develop a method to calculate the particle-hole matrix ele- 
ments of interest here. We then present the results of this analysis and make 
some conclusions in Section 3. 



2 Si/2-doublets 

Consider nuclei like O; the ground-states of such systems are, in fact, particle- 
hole states. One process by which nuclei of this type are created is the reaction 
(7r+,ir+) on target nuclei with closed proton and neutron shells [13,14,15]. 
During the course of this reaction a neutron is converted into a A. As a result, 
a neutron hole is also created which, for the ground-state, inhabits the out- 
ermost neutron shell. The angular momentum of the A and the neutron hole 
couple to form a multiplet. However, due to the fact that in the ground-state 
the A occupies the lsi/2 shell, there are only two states in these multiplets. It 
is these configurations that wc refer to as Si/2-doublets. The reaction (e, e'K^) 
is another process used to create nuclei of this type [11,12]. This process differs 
in that a proton hole is created here and that greater resolution is possible. 

In order to calculate the size of these splittings, we must first construct an 
effective interaction to model these systems. The procedure we follow here is 
similar to the method developed by Machleidt and others [19]. In this scheme, 
the effective NN interaction is represented by the exchange of mesons; then, a 
hierarchy of Feynman diagrams depicting all the possible interactions is devel- 
oped to reproduce the NN interaction. The contribution for a given diagram 
is just the product of the vertex contributions and the meson propagator, or 



where Ua are the Dirac free fields, D^q^) is the meson propagator, and the 
vertex factors, Fi and F2, are taken directly from the underlying effective 
lagrangian through the relation 



where 0/? is some meson field operator. However, nuclei are comprised of bound 
nucleons and not free fields. Therefore, in the present calculation, we improve 
on this system by replacing the free fields with the Kohn-Sham (or Hartree) 
wave functions [3] 




(1) 



5Li = iipFiipcpis 



(2) 



'4^a '4^nKmt 



r 





(3) 
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The contribution for any given diagram now takes the form (following from 
Eq. (1)) 

V{q,) = (i^J^i^a) D{q,) {i^^,T,ij^') (4) 

where we now define the effective interaction in momentum space as 

VEFTiq,.) = riD{q^)r2 (5) 

and Veft{(1ij,) — l^^f^VEFriQn)- The meson propagators relevant to this 
analysis are (here the conventions of [20] are used) 

1 1 

iql + ml 

The second term in the vector meson propagator vanishes in any calculation 
due to conservation of the baryon current (as a result, it is henceforth ne- 
glected). The relevant nucleon vertex factors are 



scalar 

vector (6) 



igs ; scalar 
-gVYn ; vector 

— ^(J^ii/g/i ; vector tensor (7) 



(Note that the A vertex factors are equivalent to their nucleon counterparts 
except for the coupling constants [6].) Now a series of diagrams can be written 
down to represent the NN interaction, each of which will contribute to this 
effective interaction in the form of Eq. (5). Once a Fourier transform is per- 
formed on the effective interaction, it can be used to calculate the two-body 
matrix elements that govern the particle-hole interactions of interest here. 

Now let us return to the case of single A- hyper nuclei. The effective AN inter- 
action is determined via the method outlined above and follows directly from 
our underlying effective lagrangian (see [6]). In this case, no single isovector 
meson coupling is allowed; as a result, we confine the following discussion to 
isoscalar, scalar and vector, exchange. As a first approximation, we include 
only the leading-order contributions arising from contact vertices; this corre- 
sponds to simple scalar and neutral vector couplings. To acquire each portion 
of the effective interaction, we simply take the product of the nucleon vertex 



4 



factor, the meson propagator, and the A vertex factor. The simple scalar and 
neutral vector exchange contributions to the effective interaction take the form 



Vvin,) = (-9KT<") (7^7^) (-9vA7f) (9) 



Next we take the Fourier transform of Eqs. (8) and (9), neglecting retardation 
in the meson propagator (i.e. ?4 — > 0), and we get 

Vsiru) - -^-^ (10) 



Tj ( . gygvA (1) (2)^:!^ 



= Vvii){ruh^M^ + ^y(2)(ri2)7?^7f ^ (H) 



where ru = {rj + r1 — 2rir2 cos ^12 )^/^. Notice that this corresponds to simple 
Yukawa couphngs to both the scalar and neutral vector mesons. The effective 
interaction to this order was used in calculations conducted in [6]; unfortu- 
nately, it proved inadequate to fully describe the ground-state splittings in 
single A-hypcrnuclci. As it turns out, the effective interaction to this order 
includes only the spin-spin force. However, this neglects the fact that the ten- 
sor force is known to play a crucial role in these systems [7,8]. Therefore, the 
natural extension of this approach is to include the higher-order contributions 
involving tensor couplings and thereby incorporate the tensor force into the 
calculation. 

There are three higher-order terms containing gradient vertices (or tensor 
couplings) of the neutral vector meson; their resulting contributions to the ef- 
fective interaction are now considered. The contributions from neutral vector 
meson exchange with a tensor coupling on one vertex to the effective interac- 
tion are 




5 



where the contribution with the tensor couphng on the nucleon (A) vertex is 
denoted by Vnt (V^at)- Notice that two terms arise in both Vnt and Vat 
due to the fact that the quantity V^^, = d^V^ — d^V^ appears in the lagrangian 
[6]. Now the Fourier transforms of Eqs. (12) and (13), neglecting retardation 
in the meson propagator, yield respectively 



tT / N ifvgvgVA my r . r / ^) 



(2) (2) 

74 75 



P^-r)7?M^^}(l + ^ 



mvri2/ ri2 

mm .jT N„,(i)„,(2) 



= Vnt^x) {n2)irir + VNTm {ri2)irir (i4) 



,2 



W TkQv my 
Att 2M 

1 \ e~^^^^'^ 



TT I X t,yTKyv inv f-r (1) / (2) Al (1) (1) 



K'•f)^^'7f}(l + — ) 



myr 12/ ri2 



= 1^AT(i) (ri2)7r'7r' + VKTi2) {ri2)lfir (15) 
where the following relation 



= X g-a(i) -571') (16) 

has been used (and also for (1) ^ (2)). It is interesting to note that in [19], the 
terms corresponding to Vnt{i) and V^Ar(i) both develop into a combination 
of spin-spin and tensor forces while the terms corresponding to V^ArT(2) and 
^AT(2) both form a combination of central and spin-orbit forces. 

For vector meson exchange with a tensor coupling on both vertices, the effec- 
tive interaction is 



T7 / ^ / fv9V (1) \ A 1 \ f9TA9v (2) 
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.(fvgv (1) \ A 1 \( 9ta9v (2) \ /T7^ 



V AM '"'^y \iql + ml)\ AM 



Taking the Fourier transform of Eq. (17), again neglecting retardation in the 
meson propagator, gives 



12n 8M2 



ri2 

= FTT(i)(ri2) (2 + 7j'^7f ) +FTT(2)(ri2) (l -T^'^T^O (!§) 

where 

5i2(f) = 3(a(^) • f) • f) - a(^) • a^^^ (19) 



Lastly, wc combine the interactions from all five contributions into a single 
effective interaction, or 



= 7? ^7? ^ (^5 + Vv + Vmt + Fat + Vtt) 

t7 (1) (2) , t7 I t7 (1) (1) (2) (2) , t7 (2) (2) 

, -FT (1) (2) , -rj (1) (1) , T7 (1) (2) 

+ VnT(2)11 74 + V^Ar(l)74 75 + ^Ar(2)74 75 

+Vrn.) (27iSf +7lSJSf7f) 

+FTT(2)(7lSf^-7lSJSf7?0 (20) 

Now that we have constructed an effective interaction, it can be used to de- 
termine the particle-hole splittings. In order to accomplish this, we must first 
calculate matrix elements of the following varieties 

((ni/iji)(n2/2j2)^M|y,(ri2)|(n3/3j3)(n4/4j4)J'M') (21) 

((niZui)(n2Z2i2) JM|y,(ri2)<T(') • B^''\n^hj^){nihji)J'M') (22) 

((ni/iji)(n2/2j2) JM|V,(ri2)(a(i) • f)\{n^h3:,){n^h3^)J'M') (23) 

{{n^kh){n2l232)JM\Vi{n2)t [(^^^^ x f) ■ a^^)] \{mhh){nihji)J' M') (24) 
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where the single-particle wave functions are specified by {nlj}, corresponding 
to either the upper or lower components in Eq. (3), and Vi(ri2) is some part of 
the effective interaction. Next, we expand each part of this effective interaction 
in terms of Legendre polynomials [5] 



Viiru) = E r2)Pk{coseu) (25) 

k=0 

= f:n{n,r2)C,{l)-C,{2) (26) 

k=0 



where 

/ 47r 

Ck, = (^^) Y,,{9, ct>) (27) 

Inverting Eq. (25) yields the expression 

2k + 1 /■ 

/fe(n,r2) = — ^ j d{cosei2)Pk{cos0,2)V,in,) (28) 



-1 



In the case of Eq. (22), the effective interaction is coupled to Pauli matrices. 
Therefore, Eq. (25) is modified to 



V^.(n2)a« • a(^) = U-^r'-'fliri,r2)xf''\l) ■ xf''\2) (29) 



Here x^xil^^ ^-^^ Ckq coupled to Pauli matrices, shown by 



X?/^ = ECfc.^i.'(fc5l5'l^lV) (30) 

qq- 



Next, we reduce Eq. (23) to the form of Eq. (21) (up to a sign). This is possible 
as the operator {a^^^ ■ f) acts only on the angular portion of the Hartree wave 
functions. Using the relation 

('^^'^•O^^U^-^^'U (31) 



the expression Eq. (23) can be rewritten in the following form 

((ni/iji)(n2/2j2) JM|\/,(ri2)(a(i) ■ r)\{n^hj^){ndAH)J'M') 
= -{{nihn){n2l232)JM\V,{ri2)\{n^[hA ^ hB\3^){nihjA)J' M') (32) 
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where liA and are the / values corresponding to the upper and lower Hartree 
spinors respectively for the ith wave function. Eq. (32) is readily generalized 
to the case (1) (2). 

Similarly, Eq. (24) can be reduced to the form of Eq. (21) (up to a factor). 
Here we employ the relation 



a 



(2) 



I(^4ji)(^2^J2)JM) 



=^/2EE i(^;^j-i)(^2k)jM)(-i)-+^-^+^ r 



X 



(^;^jiii(^(^^xf)ii/4j,)(/4j;p(2)ii(/4)j,) 



(33) 



We use [21] to further simplify the reduced matrix elements. Now we can write 
all possible matrix elements to this order in terms of Eqs. (21) and (22). 

The matrix elements in Eqs. (21) and (22) arc actually six dimensional inte- 
grals. Treating the 7-matrices as 2 x 2 block matrices operating on the upper 
and lower components of the Hartree spinors, these Dirac matrix elements, 
for each term in the interaction, are actually the sum of four separate inte- 
grals. Thankfully, angular momentum relations allow one to integrate out the 
angular dependence [21]. Eq. (21) becomes 



(21) = E(12|/^(ri,r2)|34)(-l) 

A;=0 



j2+j3+J 



J 32 jl 



' Sjj'Smm' 



x((^4)jir,(l)||(Z3^)i3)((/2^)i2||C,(2)||(^)j4) 



(34) 



and Eq. (22) becomes 



(22) = EE(12|/^(ri,r2)|34)(-l)^-^+^3+7+'=+^-A 

k=0 X 



X((/l^)jl||xf'^(l)||(/3^)j3)((/2^)j2||xr^(2)||(/4^)j4) 



J J2 jl 

1, 



> Sjj'Smm' 



(35) 



where i denotes some portion of the effective interaction. The 6-j symbols 
limit the possible allowed values of k and A. The reduced matrix elements are 
evaluated using [21] and further limit k and A. 
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Now consider the remaining two-dimensional radial integrals, where the num- 
bers are a shorthand for all the quantum numbers needed to uniquely specify 
the radial wave function [5], 

oooo 

(12|/^(ri,r2)|34) ^JJ dndr2U,{n)U2{r2)fl{n,r2)Us{n)U^{r2) (36) 





Here R{r) = U{r) / r are the appropriate radial Dirac wave functions, in terms 
of Ga{r) and Fa{r). Note that as the upper and lower Hartree spinors have 
different I values, the reduced matrix elements in Eqs. (34) and (35) must have 
the corresponding, appropriate I values. 

Using the Hartree spinor representation, the particle-hole matrix element is 
expressed as a sum of Dirac matrix elements of the types shown above [22] 

vi-M - E(2^' + 1) I '"^ ^' I {IbJ'b^W^'viamJ') (37) 
J' 3b 3i J 



No exchange term is required here, since the A and the nucleon are distinguish- 
able particles. For example, the particle-hole matrix element for the Vv{2) is 



^3'2;14(^n2)) = (-ir^^'^^'EE(-l)' 

k \ 

^^^'^ ( f dndr2 {Gi{n)F3{n)fP\n,r2)G2{r2)F4{r2) 
3i 33 J 

x((^l4)^'l|IX?''Hl)ll(^3B^)i3)((/2A^)j2||X?''H2)||(i4Bl)i4) 

-G^in)Fsin)fP\n,r2)F2ir2)G,ir2) 

x((/l4)i'l|lxf'^(l)ll(/3B^)j3)((/2B^)j2||xf'^(2)||(/44)^4) 

-F,{r,)G3{r,)fp\n,r2)G2{r2)F,{r2) 

x((^lB^)j\||xf'Hl)||(/34)i3)(a24)i2||X?''^(2)||(/45^)j4) 

+F,in)Gsin)f^^'\n,r2)F2ir2)G,{r2) 

x((/lB^)jl||X?''Hl)||(/34)j3)((/2B^)j2||xf'^(2)||(/44)j4)} (38) 

Now the sphtting, for a Si/2-doublet, is just the difference between the particle- 
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hole matrix elements of the two available states, or 

The substitutions used to get the appropriate indices for this case are n = 1, 3 
and A = 2, 4. The solution to the Kohn-Sham equations yields a single-particle 
energy level for the ground-state, E^. As previously mentioned, this level is in 
fact a doublet; however, Eq. (39) determines only the size of the sphtting. In 
order to determine the position of the doublet relative to £'a, one needs the 
relation 

^ (2 J -I- 1) 5e = (40) 
J 

We now have a framework with which to calculate the size of the Si/2-doublets 
of the single A-hypernuclei of interest here and to determine their location 
relative to The problem is reduced to Slater integrals and some algebra; 
the 6-j and 9-j symbols are determined using [23,24]. The Dirac wave functions 
needed to solve the integrals are taken as the solutions to the radial Kohn- 
Sham equations [6] . Once all the parameters in the underlying lagrangian are 
fixed, the splitting is completely determined in this approach. We also mention 
that this approach is applicable to excited states and multiplets for this class 
of nuclei. 



3 Results 

Here, we discuss the results obtained from the calculation of the ground-state 
particle-hole splittings in single A-hypernuclei by the method discussed in the 
previous section. The goal of this calculation is to evaluate 8e in Eq. (39). To 
facilitate this, it is convenient to write 5e as a sum of the contributions from 
each portion of the effective interaction, or 

5e = 5e[5] + (5e[l^(l)] + 5e[\/(2)] + be\NT{\)\ + 8e{NT{2)\ + 5e[AT(l)] 
+(5e[AT(2)] + 5e[TT(l)] + 5e[TT(2)] (41) 

where these contributions are defined in Eq. (20). As it turns out, the following 
terms cancel in the splitting 

bt\S\ = de[V{l)] = de[NT{2)] = Se[AT{2)] = (42) 
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Nucleus 


Dtate 


Levels 


V{2) 


A/r(i) 


Ar(i) 


rr(i) 


TT{2) 


X ^ 

06 


A ^ 


(lP3/2)p (1Si/2)a 


o — -1 — 


/I o 
-4zo 


- /4 


-ioo 


-1Ud8 


-ZOO 


-zUli 


16l\T 


/■T^ A/In \ 

(lPl/2)p(l'Sl/2)A 




-47d 


z8o 


-lUoz 


47d 


791 


26 




(lP3/2)p (1Si/2)a 


o — 1 — 




-57 


-14d 


-yul 


-zlz 


-looz 


16 r» 
aO 


(lPl/2)n(lSl/2)A 


n~ 1 — - 


-484 




-1U71 


485 


one; 
oU5 


zz 


A Si 


(l«5/2)n (1Si/2)a 


q+ o+ 
^G5' ^ 






-49 


ACiCi 

-490 


-15U 


1 ni 1 
-lUll 


32 Q 
A 


(2.Si/2)„(1Si/2)a 




-ZZO 


-174 


-DOl 


1 no A 


-198 


oo^in 

-zzdU 


40 no 
A 


(la3/2)n (1'5i/2)a 


1 + 9+ 


one 


Q /I 

o4 


-149 


1( ( 


ZOZ 


1 n'7 
lU/ 




(la3/2)n (lPl/2)A 




-6ib 


O 1 


97 


o o c 

-385 


on 


'71 O 

-7iz 


48 

aK 


(la3/2)p (1si/2)a 


1 + 9+ 
-^GS' ^ 


-624: 


00 


-ioL) 


111 


z47 


on 


A 


(l/7/2jn llSl/2jA 


/I "~ Q — 




-0 


1 o 
-Iz 


OOQ 
-ZZO 


-140 


-ooo 


fBb 


(1,/5/2);H1^1/2)a 




-178 


8 


-38 


187 


128 


108 


fSr 


(lfi'9/2)n^(lSl/2)A 


^G5' 


-77 








-106 


-67 


-251 


208pb 


(1^13/2)n^(l'Sl/2)A 


7g5' 6" 


-14 





-1 


-53 


-37 


-104 



Table 1 



Calculation of the sx/2-doublets (and some excited state splittings for j^^N and ^^Ca) 
in select single A-hypernuclei (in keV). Here 6e is defined by Eq. (39) and the ground- 
states are marked by GS (similarly LL denotes lower level for the excited states). 

This is true for any system in which either the A or the nucleon hole has 
j = 1/2. Thus, the total splitting for these states is 

Se = Se[V{2)] + Se[NT{l)] + 5e[Ar(l)] + Se[TT{l)] + Se[TT{2)] (43) 

It is interesting to note that these terms contribute only spin-spin and tensor 
forces; no central or spin-orbit forces survive (see Eq. (42)). The remaining 
terms can be described in the following fashion: V{2) and TT(1) are purely 
spin-spin interactions, TT(2) is purely a tensor interaction, and NT[1) and 
AT(1) are both an admixture of spin-spin and tensor interactions. Note that 
the higher-order contributions incorporate the tensor force into the calcula- 
tion; the tensor force is known to play an important role in the AN interaction 
[7,8] and does not appear in the leading-order terms. Now we determine the 
particle-hole matrix elements for each portion of the effective interaction. The 
two-dimensional integrals are calculated numerically using the Hartree spinors, 
Ga{r) and Fa{r), acquired by solving the self-consistent single-particle equa- 
tions [6]. The position of the states relative to the Kohn-Sham level is deter- 
mined from Eq. (40). All of the coupling constants used in this calculation are 
taken from [6] (specifically the sets G2, which originates from [3], and M2); 
hence, there are no remaining free parameters. 
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5 10 15 20 25 30 35 40 45 50 

2 

a 

Fig. 1. Effect of the correlation function from Eq. (44) on the total ground-state 
splitting for j^^O as a function of o? . The value of the splitting for = oo is marked. 



Nucleus 




Reference 


12-R 


< 140 


[25] 


16 Q 


+26 


[9,10] 


fSi 


< -100 


[15] 


320 

A 


< 1000 


[14] 


A Ca 


< 2200 


[13] 


208pb 


< 2000 


[15] 



Table 2 

Experimental constraints on the splittings in keV. 

The contributions from the surviving portions of the effective interaction, the 
total splitting, and the resulting level orderings for a number of ground-state 
particle-hole splittings (as well as some excited states) are shown in Table (1). 
Note that the contribution labeled V{2) was the portion of the interaction that 
was investigated in [6] ; the interaction to this order failed to reproduce either 
the correct level ordering or splitting for the ground-state of ]^0. Therefore, 
the interaction was expanded to include the higher-order gradient couplings 
described above. This expanded interaction, as given by Eq. (20), now gives 
both the proper level ordering and splitting for the ground-state of ]f O, as 
shown in Table (1). Note that the inclusion of the tensor force was crucial to 
achieve the cancellation necessary to describe the small experimental splitting 
in the ground-state of JfO [9,10], in agreement with previous work [8]. Similar 
cancellation occurs for states with ji + J 2 + vr = even (where vr denotes the 
parity of the system). Unfortunately, the splittings shown in Table (1) with 
ji + i2 + TT = odd, the ground-state of ^^S for instance, are quite large, well 
outside the known experimental error bars; the experimental constraints are 
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2 

a 



Fig. 2. Effect of the correlation function from Eq. (44) on the total ground-state 
splitting for ^^S as a function of o? . The value of the splitting for = oo is marked. 

listed in Table (2). 

The fact that the interactions take the form of Yukawa potentials here implies 
that there is some large contribution from short- distance physics that is in- 
fluencing the calculation. To correct for this problem, a correlation function 
was introduced to remove the short-distance physics from the integrals. The 
following correlation function was used here 

Vcorr = (l - e-"'<('->-^<)') (44) 

A range of was investigated for both j^^O and ^^S, the effects of which 
are shown in Figs. (1) and (2) respectively. Note that, regardless of o? ^ the 
correlation junction does not alter the level ordering of the doublet; it changes 
only the magnitude of the sphtting. Also, the cancellation that yields a small 
splitting in the ground-state of is unaffected by the correlation function. 
Thus, we can improve the splittings which were quite large while retaining the 
small splittings that resulted from cancellation. 

Technically, the proper calculation in an effective field theory such as this one 
is to choose an appropriate cutoff, then add a contact term for each portion 
of the interaction, which are essentially just constants and can be fit to ex- 
periment. However, in this case not enough data is available for nuclei in the 
range of A accessible to this type of theory. The only relevant splitting that 
has been measured is the ground-state of jfO [9,10]. Therefore, we conclude 
that at best one could claim to have a single contact term for all parts of the 
interaction. This is equivalent to a one parameter phenomenological calcula- 
tion containing a correlation function in coordinate space meant to simulate 
the proper calculation. 
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Table 3 

Calculation of the Si/2-doublets (and some excited state splittings for and j^Ca.) 
in select single A-hypernuclei (in keV) using the correlation function from Eq. (44). 
Here the value of = 0.044 was used. 

The value of the cutoff that reproduced a splitting of +26 keV in the ground- 
state of is a = 0.044; this translates into a momentum space cutoff of 
A = a^niy ~ 160 MeV. The results of calculations conducted with this cutoff 
are shown in Table (3). All of the splittings are now within the experimental 
error bars. 

Fig. (3) shows the experimental error bars from the [ti^^K^) reaction [13], 
the single-particle energy level determined from the self-consistent equations, 
the splitting and level ordering corresponding to the V(2) contribution [6], the 
splitting and level ordering determined from the expanded interaction given 
by Eq. (43) without the correlation function and with the correlation function 
where = 0.044, and the experimental doublet and level ordering [9,10] for 
]^0. Figs. (4) - (7) show, for a range of nuclei taken from Table (1), the same 
contributions as Fig. (3) except that the experimental splittings have not been 
measured for these states. One can see from Figs. (3) and (6) that the addition 
of the tensor force caused the level ordering to flip and decreased the size of the 
splitting. In addition, the correlation function has only a limited effect on the 
size of the splitting. In contrast, one can see from Figs. (4), (5), and (7) that 
the tensor force is an additive contribution to the spin-spin force, hence the 
splitting becomes large. Also, the effect of the correlation function with this 
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Fig. 3. The ground-state of is plotted here, (a) is the experimental level and 
error bars [13] along with the single-particle energy level [6], (b) is the splitting 

determined from V{2), (c) is the splitting determined from the expanded interaction 
in Eq. (43), (d) is the expanded interaction with the cutoff of = 0.044, and (e) 
is the experimental splitting [9,10]. 

particular a is to decrease the splitting size to within the known experimental 
constraints. Again, note that the level orderings for the expanded interaction 
are unaffected by the correlation function for all cases. 

The predicted ground-state of is 2"", as seen in Fig. (5). This is inconsistent 
with an analysis of the emulsion data in [26] that determined the ground- 
state spin of this nucleus to be 1. However, this theory relies on spherical 
symmetry and there is evidence that j^B is heavily deformed. Also, it should 
be mentioned that this nucleus may in fact be too small for this type of mean- 
field approach. A resolution to this discrepancy will be the subject of future 
work. 

In conclusion, wc have developed a method to calculate the doublet splittings 
of select ground-state single A-hypernuclei. This method consists of supple- 
menting the self-consistent single-particle equations by constructing an effec- 
tive interaction to simulate the residual particle-hole interaction. The form 
of the effective interaction used here follows directly from the underlying la- 
grangian. Note that this formulation of the problem contains no free parame- 
ters. Retaining only the leading-order interaction terms, this calculation was 
conducted in [6] ; this level of truncation in the residual interaction was inade- 
quate to describe either the doublet size or level ordering in the ground-state 
of ]^0. To improve on this calculation, wc included in this effective interaction 
the contributions that contained gradient couplings to the neutral vector field; 
this incorporated a tensor force into the calculation known to play a crucial 
role in these systems that did not appear at leading-order. Cancellation occurs 
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Fig. 4. The ground-state of ^^S is plotted here, (a) is the experimental level and error 
bars [14] along with the single-particle energy level [6], (b) is the splitting determined 
from V{2), (c) is the splitting determined from the expanded interaction in Eq. (43), 
and (d) is the expanded interaction with the cutoff of = 0.044. 

for the states that satisfy ji + j2 + tt = even, flipping the sign from the simple 
leading-order spin-spin interaction. However, the contributions are additive for 
the states satisfying ji+ji+T^ — odd, resulting in splittings that lie outside the 
known experimental error bars. It turns out that the integrals are dominated 
by short-distance physics; as a result, a cutoff was introduced to reduce this 
contribution. This cutoff did not effect the level orderings of any state; it did 
however, reduce the size of the splittings for states with ji + j2 + = odd to 
within the experimental constraints while simultaneously retaining the can- 
cellation the yielded small splittings in the states ji + j2 + — even. Thus, 
we obtain a realistic description of the effect of the tensor couplings on the 
doublet orderings and splittings. 
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and advice, their careful reading of the manuscript, and their helpful com- 
ments. This work is supported by Indiana University and the State of Indiana 
at the Indiana University Nuclear Theory Center, 2401 Milo Sampson Lane, 
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